Impurity induced spin-orbit coupling in graphene 
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We study the effect of impurities in inducing spin-orbit coupling in graphene. We show that the 
sp 3 distortion induced by an impurity can lead to a large increase in the spin-orbit coupling with a 
value comparable to the one found in diamond and other zinc-blende semiconductors. The spin-flip 
scattering produced by the impurity leads to spin scattering lengths of the order found in recent 
experiments. Our results indicate that the spin-orbit coupling can be controlled via the impurity 
coverage. 

PACS numbers: 81.05.Uw,71.70.Ej,71.55Ak,72.10.Fk 



Since the discovery of graphene in 2004 [l[ much has 
been written about its extraordinary charge transport 
properties 0, 0] , such as sub- micron electron mean- free 
paths, that derive from the specificity of the carbon a- 
bonds against atomic substitution by extrinsic atoms. 
However, being an open surface, it is relatively easy to hy- 
bridize the graphene's p z orbitals with impurities with di- 
rect consequences in its transport properties 0, 0] ■ This 
capability for hybridization with external atoms, such as 
hydrogen (the so-called graphane) , has been shown to be 
controllable and reversible 0] leading to new doors to 
control graphene's properties. 

Much less has been said about the spin-related trans- 
port properties such as spin relaxation, although recent 
experiments show that the spin diffusion length scales 
0,11] are much shorter than what one would expect from 
standard spin-orbit (SO) scattering mechanisms in a sp 2 
bonded system 0. fn fact, atomic SO coupling in flat 
graphene is a very weak second order process since it af- 
fects the 7T orbitals only through virtual transitions into 
the deep a bands [lj|. Nevertheless, it would be very in- 
teresting if one could enhance SO interactions because of 
the prediction of the quantum spin Hall effect in the hon- 
eycomb lattice [llj and its relation to the field of topo- 
logical insulators [121 ] . 

In this paper we argue that impurities (adatoms), such 
as hydrogen, can lead to a strong enhancement of the 
SO coupling due to the lattice distortions that they in- 
duce. In fact, it is well known that atoms that hybridize 
directly with a carbon atom induce a distortion of the 
graphene lattice from sp 2 to sp 3 [l3j . By doing that, the 
electronic energy is lowered and the path way to chemical 
reaction is enhanced. Nevertheless, it has been known for 
quite sometime 14| that in diamond, a purely sp 3 carbon 
bonded system, spin orbit coupling plays an important 
role in the band structure since it is a first order effect, of 
the order of the atomic SO interaction, A"* ps 10 meV, 
in carbon Here we show that the impurity induced 
sp 3 distortion of the flat graphene lattice lead to a signif- 
icant enhancement of the SO coupling, explaining recent 
experiments 0, @] in terms of the Elliot- Yafet mechanism 
for spin relaxation 2, T3] due to presence of unavoidable 



environmental impurities in the experiment. Moreover, 
our predictions can be checked in a controllable way in 
graphane 0] by the control of the hydrogen coverage. 

We assume that the carbon atom attached to an impu- 
rity is raised above the plane defined by its three carbon 
neighbors (see Fig. [I]). The local orbital basis at the po- 
sition of the impurity (which is assumed to be located at 
the origin, R; = o = 0) can be written as: 
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where |s), and \p x .y, z ), are the local atomic orbitals. 
Notice that this choice of orbitals interpolates between 
the sp 2 configuration, A = 0, to the sp 3 configuration, 
A = 1/2. The angle 8 between the new a orbitals and the 
direction normal to the plane is cos(#) = —A/\J A 2 + 2. 
The energy of the state | tt^ ) , e„, and the energy of the 
three degenerate states |cr a .i), e CT (a = 1,2,3), are given 
by (sec Fig. [2]): 

e„(A) = A 2 e s + (1-A 2 )e p , (2) 
e a (A) = (l-A 2 )e s /3+(2 + A 2 )e p /3, (3) 

where e s w —19.38 eV (e p « —11.07 eV) is the energy 
of the s (p) orbital [13]. At the impurity site one has 
A « 1/2 while away from the impurity A = 0. 

The Hamiltonian of the problem can be written as, 
Tt = Tin + H a + SH, where Ti. n {H a ) describes the n- 
band (cr-band) of flat graphene, and STi. describes the 
local change in the hopping energies due to the presence 
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FIG. 1: (Color online). Top: Top view of the graphene lattice 
with its orbitals. The orbitals associated with the impurity 
and lattice distortion are shown in solid black, (a) sp 3 or- 
bital at impurity position; (b) sp 2 orbital of the flat graphene 
lattice. 



of the impurity and sp 3 distortion: 

«=T,I 

a=l,2,3 

+ V™ c taO ( C <?iaQ + c a 2 aQ + C CT3 ao) + h.C.j (4) 



where 
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ci.a (cj a ) annihilates (creates) an electron at the impu- 
rity, and Cnai {ca a ai) annihilates an electron at a carbon 
site in an orbital 7r (a a ) at position R^ with spin a, e/ is 
the electron energy in the impurity, and tc-i the tunnel- 
ing energy between the carbon and impurity, Se^i^A) = 
e„{A) - e^(A = 0), and 8e a {A) = ^(A) - e a (A = 0). 
In we have not included the change in the hopping 
between a a fi orbitals (the change in energy due to the dis- 
tortion is — A 2 (e s — e p )/3) and the inter-atomic hopping 
terms. In this way, we have simplified the calculations 
and the interpretation of the results. The inclusion of 
the other terms do not modify our conclusions. 

The atomic spin orbit coupling, = A°qL • S, in- 
duces transitions between p orbitals of different spin pro- 
jection [lfj. In flat graphene (A = 0), it leads to tran- 
sitions between the n and a bands. The change in the 



ground state energy in this case is rather small and given 
by: (Af ) 2 /(^(A = 0) - e a (A = 0)) « 1(T 2 mcV 0. 
However, the perturbation described by (g]) leads to a di- 
rect local hybridization V na between the 7r and a bands 
that modifies the effective SO coupling acting on the 7r 
electrons. The propagator of n electrons from position 
Ri with spin a to Rj with spin j3 can be written as: 

(TTi,a\ (e - ft) 1 \Kj,p) ~ {lTi,a\ (e - T-Ctt)' 1 \W0,a) 

x(B k , a \Hf o \TT k , )(-K k , \ (e-n^y 1 Itt^) (6) 
where \a . a ) = [|o" 10 ) + |t7 2 o, a } + |o-3o 1 a)]/v / 3anda J - iQ ,) = 
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\^3j,a)]/V^ where 4> = 27r/3. 
The propagator in ([6]) can be understood as arising from 
an effective non-local SO coupling within the 7r band 
which goes as: 



A^(0,i) « ^(ao, Q Ke-W ff )" 1 |^, Q )A° 



(7) 



which allows us to estimate the local value of the SO 
coupling as: 



A^Z{l-A 2 ). 



(8) 



As shown in Fig. [2] the value of the SO coupling depends 
on the angle (i.e., the value of A) associated with the 
distortion of the carbon atom away from the graphene 
plane. Notice that for the sp 2 case (A — 0) this term 
vanishes indicating that SO only contributes in second 
order in A"*, while for the sp 3 case (A = 1/2), the SO 
coupling is approximately 75% of the atomic value (w 7 
meV). Also observe that the dependence on the distance 
from the location of the hydrogen atom is determined by 
the Green's function G (J (0,Rj)) = (ao|(e — Hcr) _1 |CTj). 
This function, calculated for the simplified model of the a 
bands discussed in ref. (lo| , shows a significant dispersion 
in Fourier space, ranging from a maximum at the T point 
to zero at the K and K' points. Hence, the range of 
G(j (0, R) should be of the order of a few lattice constants. 

Based on the previous results we can now calculate 
the effect of the impurity induced SO coupling in the 
transport properties. Firstly, we linearize the ir band 
around the K and K' points in the Brillouin zone and 
find the 2D Dirac spectrum Q: ej-.k = ±vpk where vf 
(w 10 6 m/s) is the Fermi-Dirac velocity. In this long 
wavelength limit the impurity potential induced by ([7|) 
has cylindrical symmetry and we can use a decomposition 
of the wavefunction in terms of radial harmonics 119. 12(1. 
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22j, |23j. A similar analysis, for a system with SO 
interaction in the bulk has been studied in ref. Q. We 
describe the potential scattering by boundary conditions 
such as one of the components of the spinor vanishes at 
a distance r = R\ (of the order of the Bohr radius) of 



the impurity 24|. A Rashba-like SO interaction exists in 



3 



the region R\ < r < R 2 (region I), and there is neither 
potential nor spin orbit interaction for r > R 2 , region II 
(R 2 h of the order of the carbon-carbon distance). 

The wavefunctions in region I can be written as a su- 
perposition of angular harmonics: 
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where | ]) and | J.) are the spin states. The functions 



J n (x), Y n (x) are Bessel functions of order n, and: 



(10) 



c± = ] /l/2±AU(^v 2 F kl + (AU2r) (11) 



4 = ^l/2±AU(^v 2 F k*_ + (AUm (12) 

e is the energy of the scattered electron (k± is defined 
through ([TO])). 

The wavefunctions outside the region affected by the 
impurity, r > R 2 , can be written as: 



( Y n (kr)e mS \ |\ 
+ L ^\iY n+l {kr)e i ^+ 1 ) 6 ) 1/ 

+ Ll \iY n+2 {kr)e^+^ ) 1/ 
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and: e = vpk. The boundary conditions at r = Ri and 
r = R 2 lead to the equations: 



c + A+ J n (k+Ri) + c+B + Y n {k + Ri) + c'_A_J n (k-Ri) + c'_B_Y n (k-Ri) = 
c-A + J n+1 {k+Ri) +c-B+Y n + 1 (k+R 1 ) + c' + A-J n+1 {k-Ri) + d + B-Y n+1 {k-R{) = 

c+A+ J n (k+R 2 ) + c+B + Y n (k+R 2 ) + d_A_J n {k-R 2 ) + c'_B_Y n (k-R 2 ) = J n {kR 2 ) + C^Y n (kR 2 ) 
C-A + J n+1 {k+R 2 ) + c-B + Y n+1 {k+R 2 ) + c' + A_J n+1 {k-R 2 ) + c' + B_Y n+1 (k-R 2 ) = J n +i{kR 2 ) + C T F n+1 (A:i? 2 ) 
C.A+ J n+1 (k+R 2 ) + c-B + Y n+1 (k+R 2 ) - c' + A-J n+1 (k-R 2 ) - c' + B-Y n+1 (k_R 2 ) = C^+^kR^ 
c+A + J n+2 (k + R 2 ) + c+B+Y n+2 (k + R 2 ) - c'_A_ J„ +2 (fc_i? 2 ) - c'_B_Y n+2 {k.R 2 ) = C^+^kR^ (14) 



These six equations allow us to obtain the coefficients 
A±, B±,Cf and Cj_. In the absence of the SO inter- 
action, we have ^4+ = A-,B+ = £>_,Cj = and 
C T = -J n {kR 1 )/Y n {kR 1 ). 

We show in Fig. [3] the results for the cross section for 
spin flip processes, determined by \C\\ 2 /kF- The main 
contribution arises from the n — channel. For compari- 
son, the elastic cross section, calculated in the same way, 
is a e i w kp 1 . This is about three of magnitude larger 
than the spin-flip cross section due to the spin orbit cou- 
pling. Hence, the spin relaxation length is 10 3 times the 
elastic mean free path We obtain a mean free path 
of about 1 /xm, in reasonable agreement with the exper- 
imental results in ref. Q- This value depends quadrat- 
ically on Al (A). For a finite, but small, concentration 



of impurities, our results scale with the impurity concen- 
tration and hence the spin flip processes should increase 
roughly linearly with impurity coverage in transport ex- 
periments in systems like graphane [6| . 

In summary, we have shown that the impurity induced, 
lattice driven, SO coupling in graphcne can be of the 
order of the atomic spin orbit coupling and compara- 
ble to what is found in diamond and zinc-blend semi- 
conductors. The value of the SO coupling depends on 
how much the carbon atom which is hybridized with 
the impurity displaces from the plane inducing a sp 3 hy- 
bridization. We have calculated the spin-flip cross section 
due to SO coupling for the impurity and shown that it 
agrees with recent experiments. This results indicates 
that there are substantial amounts of hybridized impu- 
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FIG. 2: (Color online) .Top: Energy (in eV) of the tt (blue) 
and a (red) bands as a function of A according to Q; Bottom: 
Relative value of the SO coupling at the impurity site relative 
to the atomic value in carbon as a function of A according to 
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FIG. 3: (Color online). Cross section for a spin flip process 
for a defect as described in the text. The parameters used are 
Ri = lA R 2 = 2k and A' ao = lmeV (blue) and Af D = 2meV 
(red) . 



rities in graphene, even under ultra-clean high vacuum 



conditions. Experiments where the impurity coverage is 
well controlled can provide a "smoking-gun" test of our 
predictions. 
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